Expansion of Solid Helium into Vacuum: the Geyser Effect 

Robert E. Grisenti^'^, Manuel Kasz^, J. Peter 
Toennies^, Giorgio Benedek^, and Franco Dalfovo^ 
^ Max-Planck-Institut fiir Stromungsforschung, 
Bunsenstr. 10, 31013 Gottingen, Germany 
"^Institut fiir Kernphysik der J. W. Goethe Universitdt, 
August- Euler-Str. 6, 6O486 Frankfurt /Main, Germany 
^ Dipartimento di Scienza dei Materiali and INFM, 
Universitd di Milano-Bicocca, Via Gozzi 53, 20125 Milano, Italy and 
^ Dipartimento di Matematica e Fisica and INFM, 
Universitd Gattolica, Via Musei 4h 
25121 Brescia, and GRS-BEG, Trento, Italy 

Abstract 

The particle flux through a two micron diameter orifice into vacuum from a source chamber fihed 
with sohd ^He exhibits a striking periodic behavior similar to that of a geyser. This phenomenon 
is attributed to a periodic collapse of the solid inside the source induced by the accumulation of 
excess vacancies injected at the orifice. The fiux-time curves agree well with a kinetic model and 
provide direct information on the diffusivity of vacancies in solid ^He. 

PACS numbers: 64.80.-s, 67.80.Mg 
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Lattice vacancies are always present in an otherwise perfect crystal in thermal equilibrium 
because of the corresponding increase in entropy. Regarded classically, these vacant sites 
are localized defects which move only occasionally in a thermally activated diffusion process. 
In solid helium, however, as a consequence of the light atomic mass and weak attractive 
interaction potential the particles execute zero-point oscillations which cover a large fraction 
of the cell sizeQ]. Owing to the large amplitude of the zero-point motion tunneling of a 
vacancy into a neighbor lattice site is expected to occur at significant rates. Vacancies thus 
become delocalized entities which can be treated as elementary quasi-particle excitations of 
the crystal occupying a finite band of energy states 3, 0, 1^ , which is a unique feature of 
quantum solids. 

Despite several decades of intensive experimental effort the present understanding of 
the underlying mechanism of vacancy diffusion in quantum solids is still incomplete, as 
demonstrated by the variety of conflicting data for the basic vacancy parameters (see P, 0] 
and references therein). Very likely this arises from the difficulties in the interpretation 
of results obtained by indirect methods, especially in the absence of a realistic model of a 
vacancy^l,!^. Furthermore, the previous experiments afforded little control over the vacancy 
concentration, which drops exponentially with decreasing temperature and is thus probably 
too low to reveal significant quantum effects. 

This letter reports on a quite surprising discovery related to the above issues. In the 
expansion of solid ^He into a vacuum system through a small pinhole orifice the particle flux 
exhibits a periodic sequence of Geyser-like bursts. The periods are in the range of seconds 
or minutes, depending on the pressure and temperature. The phenomenon is interpreted in 
terms of excess vacancies injected into the solid at the orifice which diffuse inwards eventually 
causing a periodic collapse of a macroscopic portion of the solid, thereby resetting the initial 
conditions. The observed temporal behavior is shown below to be closely related to the high 
mobility of the vacancies, and provides direct information on vacancy diffusion in solid ^He. 

A standard molecular beam source in which the gas is replaced by solid ^He is used in 
these experiments. This simple arrangement represents an idealization of previous studies of 
plastic flow of solid ^He for the determination of vacancy difrusivity0,0. In the apparatus, 
shown in Fig. 1, the helium initially in the gas state is admitted into the source at a pressure 
Pq and a temperature Tq maintained by a liquid- He bath cryostatjl^. The temperature was 
determined from the He vapor pressure above the liquid helium bath. The pinhole orifice 
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with cross section Sex has a diameter d = 2 /im and a comparable length, and is located 
coaxially at the end of a long 1 mm-dia. copper tube with cross section 5*0. After passing 
through a separately pumped source chamber the narrow beam reaches the detector chamber 
through a 2.5 cm-dia. aperture, located about 10 mm beyond the orifice. The beam flux is 
proportional to the detector chamber pumping speed times the measured pressure increase 
APd in the detector chamber, which is between ~ 3.5 • 10~^ and ~ 6 • 10~^ mbar. 

In experiments at a constant source temperature Tq [Fig. 2(a)], the beam flux is measured 
as the source pressure Pq is gradually increased stepwise from below to above the melting 
pressure Pm(Tn)- At Pq < -^111(^0) the liquid flows at a constant rate steadily increasing with 



increasing Pq 



As soon as Pq exceeds Pml^o) the flux exhibits bursts with period tq, each 



burst being characterized by a sudden increase in intensity followed by a slowly decaying 
tail. The period tq is found to increase with increasing source pressure Pq [see inset of Fig. 
2(a)] and, as shown in Fig. 2(b) at a given temperature and pressure, is constant within 
< 5%, even in experiments over several hours. 

The plots of the period tq as a function of pressure at constant temperature [inset in 
Fig. 2(a)] or of temperature at constant pressure (Fig. 3) show that the period vanishes at 
either the melting pressure Pm or the melting temperature T^. The critical points at which 



To vanishes follow, within the experimental error, the me 



ting curve (see inset of Fig. 3). 



Optical observations like those reported in Ref. [llj indicate clearly a narrow liquid 
jet at the nozzle exit for all of the source conditions considered here. This indicates that 
there must be a solid-liquid (s/1) interface with an effective cross section Ssi (< 5*0) at 
some intermediate distance Tsi from the nozzle where the pressure approaches the melting 
pressure Pm, as shown schematically in Fig. 4. In the liquid layer r^i < r < the pressure 



obeys Bernoulli's equation 



111 ], with the average velocity of the liquid increasing from the s/1 



interface to the orifice exit inversely proportional to the reduction in cross section 5'ex/5'si. 
The exit liquid velocity is simply related to APd by Wex = CAPDVa^t/cP ~ 10^ m/s, where 
Vat is the atomic volume of the liquid and C is an apparatus constant. 

The near equality of the time averaged fluxes above [broken line in Fig. 2(a)] and below 
the melting point suggests a plastic flow of the solid inside the source cylinder similar to 
that of a liquid jl2|. A crucial difference between the dynamics of the solid and liquid, how- 
ever, is expected in the region where the flux cross section shrinks from the source cylinder 
area 5*0 to ^si. In this region the plastic flow of the solid is ensured by two mechanisms 
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involving the motion of defects: one is the plastic deformation allowed by the motion of dis- 
locations; the other is the injection into the solid of vacancies generated at the s/1 interface 
with a corresponding release of atoms from the solid into the liquid. Driven by their inward 
relaxation [l3| the vacancies climb up the strong pressure gradient. Although many vacan- 
cies segregate at extended defects, e.g., at the axis of edge dislocations causing dislocation 
climbing a large fraction will be pinned in the upstream region at tq where the pressure 
gradient drops to zero and the pressure quickly approaches Pq (see Fig. 4). While disloca- 
tion climbing will support a constant flow, the accumulation of excess vacancies around tq 
is held responsible for the periodic bursts in the beam flux intensity. A periodic behavior of 
the particle flux due to a stick-slip motion of the solid inside the source is ruled out since 



in this case one would expect a decrease of tq with increasing Polla]; whereas in the present 
experiments the period increases with increasing source pressure [see inset of Fig. 2(a)]. 

To model the oscillatory part of the flow it is assumed that a percolative collapse, occur- 
ring at some critical vacancy concentration X^, leads to a bleaching of vacancies and resets 
the initial conditions. By equating the total flux of the liquid at the orifice to that at the 
s/1 interface, the time dependent liquid exit velocity can be written as 

Ue.it) = iS,l/Se.) [Uo + KtJv(t)] , (l) 

where uq is the stationary velocity due to the plastic flow of the solid inside the source tube 
and the second term accounts for the time dependent vacancy current density jv{t) at the 
s/1 interface, jv(^) consists of a diffusive and a drift term 

jv(t) = -D^x'^{rsi, t) + u^x^{r^i, t) , (2) 

where Xv{r, t) is the excess in the vacancy concentration with respect to the local equilibrium 
value in the thin solid layer of thickness i = vq — Vsi (Fig. 4) and x^r, t) its space derivative. 
Xv{r,t) is determined by solving the kinetic equation in one dimension jiol 

(9Xv „ d Xy dx^ Xy , r~i f j.\ /o\ 

^ = ^v^-«v^-- + G(r,t). (3) 

In Eqs. (2) and (3) the vacancy diffusion coefficient is related to the mobility /iv through 
the Nernst-Einstein relation = kTfi^, and = fJ^vF — uq is the vacancy drift velocity 
with respect to the laboratory frame under the action of a force F. The inward relaxation 
of the lattice surrounding a vacancy is characterized by a local change AV^ < in the 



vacancy volume and causes a net force field F ^ —AV^VP on the vacancy driving it 
towards higher pressures. Thermal generation and recombination processes are included 
in Eq. (3) through the vacancy recombination time Tj.. The external generation function 
G{r,t) = X^[—6(t)9{r — rsi)9{rQ — r) +Us6{r — rsi)9(t)], where 9 is the Heaviside step function, 
consists of two terms. The first term models the sudden depletion of vacancies by an amount 

occurring in the narrow region < r < tq at time t = 0. By assuming a depletion of all 
the vacancies, X° is approximately given by an average equilibrium vacancy concentration 
of ~ 10^° cm"^!^!. The second term describes the subsequent vacancy injection at the 
s/1 interface with a vacancy generation velocity UsH]. 

Equations (2) and (3) predict in good agreement with the experiment the temporal be- 
havior of the flux within a single period. For simplicity, in Eq. (2) has been assumed 
to be constant in the solid layer rgi < r < vq. Figure 2(b) compares the calculated and 
experimental flux curves for Tq = 1.74 K and Pq = 31 bar. The best fit was obtained with 

— 1-3 • 10~^ cm^/s f corresponding to /iv ^ 5.4 • 10^° s/g), Mv ^ 2 • 10^^ cm/s, Tj. 60 
s, and with = 2mv|12|. The order of magnitude of D^, which is consistent with earlier 
measurements is typical of a liquid, thereby indicating non-classical diffusion of vacancies 
in the solid. 

Since the liquid is incompressible the exit velocity Uex is related to the pressure of the 
liquid at the s/1 interface, P^i, by the BernouUi equation ~ (2Vat-Psi)^''^/"^5 where m is 
the mass of a He atom. Thus, the oscillations in the detected signal APd correspond directly 
to small oscillations of the pressure at the s/1 interface Psi [Fig. 2(a)], right-hand ordinate 
scale]. The oscillations of Psi with respect to Pm are in turn a natural consequence of the 
variable concentration of vacancies in the solid which affects the local melting pressure. 

The excess critical vacancy concentration = X^ — X^ causing the collapse of the solid 
can be estimated from the measured Tq. Within the present kinetic model the period tq is 
approximately given by tq ~ r^{xl/ X^)(£ / icY^'^ , where = AD^/ul and = iDv^rY^^ is 
the diffusion length. By assuming AVv ~ —60 A^|13], it follows that £ fa 40 fim and hence 
xl/X^ 10. Since the integrated intensity over a full period is a measure of the amount of 
injected vacancies, the observed growth of the oscillation amplitude with increasing source 
pressure [Fig. 2(a)] suggests that at higher Pq a larger excess vacancies concentration is 
needed to induce the collapse of the solid in the source. This is consistent with the expected 
decrease of the equilibrium vacancy concentration with increasing pressure j^. 



5 



In conclusion, the temporal behavior of the beam intensity in expansions of solid ^He 
into vacuum through a pinhole orifice provides a remarkably direct way to explore vacancy 
diffusion and mobility. The orifice acts as a well defined source of vacancies, which gradually 
diffuse into the solid and accumulate in the region where the local pressure approaches the 
externally applied source pressure. Periodically, the concentration of vacancies in this region 
reaches some critical value at which the solid collapses, resetting the initial conditions. The 
extreme simplicity allows this type of experiment to be easily extended to a systematic 
study of vacancy diffusion in systems other than helium. A particularly intriguing case is 
represented by solid p-B.2, for which previous studies suggest at T 10 K the transition 
from a thermally activated to a purely quantum mechanical vacancy diffusion process 3] ■ 
The present experiments are also ideally suited to address the exciting possibility of Bose- 
Einstein condensation of vacancies in bulk solid ^Hei 



19 



20|. Despite many attempts [21|, 



22l |. perhaps because. 



only indirect evidence has so far been found for this "supersolid" state 
if it exists, the equilibrium vacancy concentration is too small to make it detectable. In the 
present experiments a much larger vacancy concentration is achieved under non-equilibrium 
conditions, thereby enhancing the possibility of BEC in solid ^He. Indeed, the probable 
observation of a supersolid phase in solid helium confined in Vycor, a configuration that 
is likely to be more heavily populated with vacancies than bulk helium, has recently been 
reported 2^ 

Useful discussions with L. Reatto, F. Pederiva, and W. Schroter are gratefully acknowl- 
edged. We thank R. Dorner for carefully reading the manuscript. R. E. G. and M. K. thank 
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FIG. 1: Schematic diagram of the source and vacuum system. Helium gas is admitted at a 
pressure Pq into the source tube inside a copper block which is attached to a liquid-He cryostat 
at a temperature Tq. The liquid jet passes into the source chamber at a pressure of ~ 10~^ mbar 
into the detector chamber. The inset shows a magnification of the source region. The 2 /xm-dia. 
pinhole is centrally located in a thinned region of an 125 fj,m- thick 2 mm-dia. nickel platelet (4) 
at the end of a 1 mm inner dia. copper tube (5). The gold ring (3) serves as a pressure seal. The 
stainless steel retainer (2) is compressed by the outer stainless threaded cap (1). 
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FIG. 2: (a) The detector chamber pressure increase APj) as a function of time t for Tq = 1.87 K. 
The source pressure Pq is gradually increased stepwise from below to above the melting pressure 
-fm(^o) ^-s indicated by the marks in the upper abscissa scale. The right-hand ordinate scale 
indicates the calculated pressure of the liquid at the solid/liquid interface corresponding to the 
observed beam intensity. The inset shows explicitly the oscillation period tq as a function of Pq. 
(b) Four complete periods measured at Pq = 31 bar and Tq = 1.7 A K are enlarged to show the 
recurring invariant shape with a constant period (tq = 82 s). The period and shape are very well 
reproduced by the pulse shapes, shown below the experiments, predicted by the excess vacancy 
transport theory Eq. (3). 
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FIG. 3: The measured oscillation period tq as a function of source temperature Tq at Pq = 32 bar 
(circle), 35 bar (square), and 41 bar (diamond). In the inset the points in the P — T plane where 
the observed period tq vanishes all lie on the melting curve of the ^He phase diagram. 
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FIG. 4: (a) Schematic diagram showing the flow tubes near the orifice inlet and (b) the dependence 
of the pressure as a function of the distance r from the orifice for the liquid (r < rgi) and for the 
solid in equilibrium with the liquid at the beginning and end of each pulse indicated by the cross- 
hatched region. In the solid layer rs\ < r < rg the flux lines start converging towards the orifice 
producing a strong pressure gradient. 
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